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ABSTRACT 



We analyze the structure of higher-order radiative corrections for processes with unstable par- 
ticles. By subsequently integrating out the various scales that are induced by the presence of 
unstable particles we obtain a hierarchy of effective field theories. In the effective field theory 
framework the separation of physically different effects is achieved naturally. In particular, we 
automatically obtain a separation of factorizable and non-factorizable corrections to all orders 
in perturbation theory. At one loop this treatment is equivalent to the double-pole approxi- 
mation (DPA) but generalizes to higher orders and, at least in principle, to beyond the DPA. 
It is known that one-loop non-factorizable corrections to invariant mass distributions are sup- 
pressed at high energy. We study the mechanism of this suppression and obtain estimates of 
higher-order non-factorizable corrections at high energy. 



1 Introduction 



An important area in the physics programme of a future electron-positron hnear colhder (LC) 
belongs to the detailed studies of the pair production of heavy unstable particles (W and Z 
bosons, top quarks, SUSY particles, etc.) One of the main goals of these measurements is 
the high precision determination of their parameters, primarily their masses. 

In recent years there has been a growing interest in high energy photon colliders (see e.g. |^), 
which in some cases offer certain advantages over the e~^e~ collisions in the exploration of the 
unstable particles, see for example . Pairs of unstable particles can also be produced in hadron 
collisions at the Tevatron and LHC All these studies would allow a variety of detailed tests 
of the Standard Model to be performed, both in the electroweak and in the strong interaction 
sectors. This, in turn, requires very accurate theoretical knowledge of the production and decay 
mechanism and of interference effects between them. This interference is caused by the large 
width, r, of many of these objects. The widths of the W, Z and t are all of the order of 
2 GeV, and many supersymmetric and other beyond-the-Standard-Model particles would also 
have widths in the GeV range. 

Although many results of this paper are applicable to the non-relativistic case as well, we 
concentrate on relativistic and ultra-relativistic energy regimes. In the center-of-mass frame 
the physical picture of processes involving unstable particles can be viewed as a consequence 
of several subprocesses. First the unstable particles are produced during a short time scale 
0{1/E). Then during a large proper time, ~ l/F, they separate from each other on a large 
distance, ~ l/F ■ E/M. Finally, they decay during a short time, ~ l/E, into their decay 
products. Here, E, M and F are the energy, mass and width of the unstable particle respectively. 
If E^ 3> MF, the production and decay subprocesses are well separated from each other. The 
so-called pole-scheme |^ is a framework which makes use of the presence of two scales in the 
problem by expanding the amplitudes in F/M. In particular, one can keep only the leading 
terms, neglecting all the terms suppressed by powers of F/M. In the case of pair production 
of unstable particles such an approximation is usually called the double-pole approximation 
(DPA). The DPA guarantees gauge invariance of the calculation and simplifies it significantly. 
The DPA allowed the calculation of the full 0{a) electroweak corrections to the pair production 
of W and Z bosons in e^e~ collisions 0. It is essential for the DPA that unstable particles 
are close to resonance and that E"^ ^ MF. This is also the regime which we will specialize to 
throughout this paper. 

Within the DPA at one loop one can readily classify the radiative effects into two types: 
factorizable, which act inside the separate hard subprocesses (production and decay); and non- 
factorizable, which interconnect various hard subprocesses. Depending on the subprocesses they 
interconnect, the non-factorizable corrections can be of two types: decay-decay and production- 
decay. The non-factorizable corrections are believed to to be an important ingredient of the 
theoretical predictions. In particular, they affect the measurement of the mass of the unstable 



particle. They have been a subject of intensive study recently, see for example 0, [|6|-pl| 



The properties of perturbative non-factorizable corrections depend crucially on the level of 
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"inclusiveness" of the distributions. The distributions can be of three distinct types. For vir- 
tual corrections only two of them are relevant: invariant mass and angular distributions. Since 
real gauge-boson (photon, gluon) radiation contributes to the non-factorizable corrections as 
well one can consider additional distributions where the non-factorizable real radiation is added 
and an integration over the real emission phase space is performed. At one loop in QED the 
non-factorizable corrections to the distributions inclusive in real photon and invariant mass are 
suppressed by at least aF/M, [0, In [|10| general arguments are given indicating that 



this should still be the case at higher loops. As a consequence, for a particular class of the 
one-loop non-factorizable corrections (initial-final state interferences) the virtual contribution is 



cancelled by the real emission even for the distributions exclusive in the invariant mass. In [12 



a first calculation of one-loop non-factorizable corrections exclusive in invariant mass but inclu- 



sive in real radiation was performed. In |13, 14 1 the one-loop non-factorizable corrections were 



calculated for completely exclusive distributions. Although the latter calculations were origi- 
nally performed for pair production of W bosons, the results were applied also to other cases, 
for example to one-loop QCD non-factorizable corrections in top-quark pair production |T^ 
and to one-loop QED non-factorizable corrections to the pair production of Z bosons 



Even though one-loop non-factorizable corrections are extensively studied by now, almost 
nothing is known about the higher-loop non-factorizable effects. An explicit calculation of two- 
loop non-factorizable corrections is well beyond present capabilities. Even a precise definition 
of higher-loop non-factorizable corrections is lacking at present. Higher-order non-factorizable 
corrections can be quite important, especially in QCD. Note that all of the one-loop results 
above did not assume the distributions to be inclusive in the angles. In [1^ it was found 



that the one-loop non-factorizable correction to the distribution inclusive in real photons and 
angles, but not in the invariant mass, is suppressed by a(M/£')^ at high energies, E 3> M. 
This is consistent with an earlier observation of Ref . , where a similar effect was observed 
for a simplified toy model. A similar energy dependence is also found for the one-loop QCD 



interconnection effects in tt production |jT5|, |T3- It is the main purpose of this paper to study 



the mechanism of the high energy suppression. This allows us to obtain information about the 
higher-order non-factorizable corrections to the distributions inclusive in angles. 

An indication of the importance of higher-loop QCD corrections comes from the studies 
of non-perturbative QCD interconnection effects. These effects are an essential source of the 
systematic uncertainties in the reconstruction of the H^-boson or top-quark parameters. Recall 
that not far from threshold the typical decay time, r ~ 1/F O.lfm, is much shorter than the 
characteristic hadronization time, rhad ^ 1 fm. Thus, the final state hadronic systems overlap 
between pairs of resonances (VF^VF~, Z^Z^, ti, etc). Non-perturbative effects at the hadroniza- 
tion stage are usually modelled by a colour rearrangement between the partons produced in the 



two resonance decays and the subsequent parton showers [|1^. This topic is deeply related to 
the confinement physics. Thus, interconnection studies not only offer an opportunity to inves- 
tigate the dynamics of unstable particles, but they also open new ways to probe confinement 
forces in space and time |T8|, |19[. 

So far, the experimental analysis of the non-factorizable effects have been performed only 



at LEP2 and mainly in the context of precise W mass measurements [2C]. A LC will allow a 
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series of very accurate measurements of W~^W~ , ZZ, and tt production in a wide energy range. 
Especially promising for the purposes of precise mass determination and interconnection studies 
look high luminosity runs of a LC in the threshold regions ||19|, |2l|] . A detailed knowledge of the 
energy dependence of the interconnection effects in W^W~ production would allow to choose 
the optimal strategy for their studies at a LC. 

The purpose of this paper is twofold. First, we discuss the separation between the different 
radiative phenomena to all orders in perturbation theory. We find the following different effects: 

• There are factorizable and non-factorizable corrections. 

• The non-factorizable corrections can be of the production-decay and decay-decay types. 

• The non-factorizable corrections of each type can receive corrections due to two effects: 
interaction between the production/decay dipoles and propagation corrections. 

We start with the separation of these effects in QED within the DPA by analyzing Feynman 
diagrams. This is an extension of the well studied one-loop QED case. Subsequently we find 
an effective field theory interpretation of this separation, based on the presence of a hierarchy 
of scales in the problem, E"^ ^ MT. The existence of an effective field theory allows us to 
generalize the separation between different effects to the more complicated (non-abelian) QCD 
case. Also the DPA appears naturally within the effective field theory framework. This implies 
- at least in principle - the possibility to study F/M suppressed (beyond the DPA) effects in 
a consistent way. 

Second, we focus on the high-energy, E ^ M, behavior of the non-factorizable correc- 
tions to the distributions inclusive in angles of the decay products. We study the mechanism 
responsible for the suppression of the non-factorizable corrections at high energy. The sup- 
pression mechanism works differently for the different effects we identified. In this paper we 
limit ourselves to the decay-decay interferences, leaving aside the more complicated case of 
the production-decay interferences. We illustrate the suppression mechanism by estimating the 
QED decay-decay interferences, using a high energy expansion of Feynman diagrams. Subse- 
quently we find an effective field theory interpretation of this mechanism. Based on the effective 
field theory approach we develop a general framework that allows us to generalize the estimates 
to the non-abelian QCD case. Finally, we consider some topical applications of our results. 

The paper is organized as follows. In Sect. ^ we discuss the separation of different effects 
in QED by analysing Feynman diagrams. We also derive the high energy estimates of the 
non-factorizable decay-decay interferences. In Sect. |^ we analyse the structure of the modes 
induced by the hierarchy of scales, E and F. By subsequently integrating out different modes 
we construct an effective field theory that provides us with a separation between the various 
radiative effects together with the high energy estimates. In Sect. ^ we discuss some physical 
applications of the high energy estimates of the non-factorizable corrections. Sect. ^ gives our 
conclusions and outlines the possible future developments and applications of our results. 
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2 Estimates from high energy expansion of amplitudes 



In this section we investigate higher-loop QED non-factorizable corrections in a toy model by 
considering the high-energy expansion of Feynman diagrams. We work within the DPA, neglect- 
ing all contributions suppressed by F/M. The separation of factorizable and non-factorizable 
corrections is achieved by using the non-factorizable currents. This procedure is an exten- 
sion of the one- loop separation procedure p, |14|. We subsequently study non-factorizable 
corrections to the distributions inclusive in decay angles in the high energy limit, E ^ M. 

2.1 Separation of various contributions 

In the Born approximation the hierarchy of scales E"^ ^ MF implies the factorization of the full 
matrix element into a product of matrix elements corresponding to hard (production and decay 
of unstable particles) and soft (propagation of unstable particles) sub-stages of the process. 
The factorization holds with an accuracy of at least (9 (F/M). 

Radiative corrections are due to radiation of photons/gluons with typical energy Q. If the 
radiation is hard, Q ^ E, the inverse propagation distance of the unstable particles, ~ T -M/E, 
remains to be the only soft scale in the problem. As a result, the factorization for the process 
works in the same way as in the Born approximation. On the other hand, if the radiation is 
soft, f2 ~ F ■ M/E, there are two soft scales. The hard matrix elements still factorize but the 
propagation of the unstable particles and the interactions via soft particles can mix now. In 
general, a part of the corrections will have a form where the propagation subprocess is factorized 
in the same way as it is factorized in the Born approximation, and a part of the corrections 
will not have such a factorized form. 

A correction is factorizable (non-factorizable) if the propagation sub-stage factorizes (does 
not factorize)with respect to the rest of the process in the same way as in the Born approx- 
imation. Corrections due to the exchange of hard particles are always factorizable whereas 
corrections due to the exchange of soft particles can be both factorizable and non-factorizable. 
In the former case the split up between factorizable and non-factorizable corrections is unam- 
biguously defined on the basis of the energy of the exchange particles. In the latter case the split 
up between factorizable and non-factorizable corrections can only be done by comparing cor- 
responding matrix elements to that in the Born approximation. One expects that factorizable 
and non-factorizable contributions factorize with respect to each other beyond one loop. 

The toy model we will consider has a neutral, scalar, unstable particle field with mass M 
and width F. This particle decays into massless, charged fermions, ip. The charged fermions 
couple to a U{1) gauge field, A^. The Lagrangian of the model is 

^ = \'P{p^ - M')0 + \f^uF^'' + vxcjxp + (i)^vij. (1) 

Here p is the momentum operator and p is the covariant momentum operator. The first three 
terms in Eq. (|l]) describe the propagation of 0, ip, and fields. The V term describes the 
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production of a pair of unstable particles from some neutral source X, whereas the V term 
describes their decay into the fermions. We will assume that all these terms can be treated as 
a perturbation. To see how factorization occurs let us consider the process 

X ^ ^{pMp2) i;^{h)i;-{k[)i;-{h)i;^{k',). (2) 

where Pi = ki + k[, i G {1, 2} and 

pI-M'^ ^ TM, pI-M'^ ^ TM (3) 

We assume there are no charged particles in the initial state. 



2.1.1 One-loop order 



The complete one-loop correction consists of one-loop factorizable and one-loop non-factorizable 
corrections. The separation between these two contributions is known . In order to introduce 
notations let us recall briefly how it works. The hard modes contribute exclusively to the 
factorizable corrections, thus the separation needs to be performed only for soft photons. Let 
us start with the non-factorizable corrections. Since charged particles appear only in the final 
state, as decay products there will be only decay-decay interference. The decay-decay non- 
factorizable correction is given by 







+ 



where we used a short hand notation for the loop integral 




(4) 



(27r)4/2 



(5) 



In Eq. A4q is the Born matrix element in the DPA. The gauge invariant currents J'12 
called the non-factorizable currents. They are given by 



J2' 



Jl\l,D,) = +e 
J^{l,D2) = +e 



Jh~Jk{ 
^2 



22 



Di + 2lpi 
D2 



(6) 



D2 - 2lp2 



where 



Di,2 =pl2-M^ + iVM, 



(7) 

is the off-shellness of the unstable particle. By inspection of the integral Eq. @ it can be seen 
that the energy of the photon relevant for non-factorizable corrections is indeed Q ~ TM/E. 
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The contribution of a hard photon is suppressed by 0{M'^T'^ / E"^) with respect to the leading 
contribution. Thus, the corresponding contributions are beyond the DPA. 

The factorizable correction is given by 

; + • + ... = All. (8) 

The important property oi M.i is that it has exactly the same dependence on Di^2 as the 
Born matrix element M.^. This is the crucial difference between the factorizable and non- 
factorizable corrections. The non-factorizable corrections contain an additional strong D12 
dependence. Thus, contrary to the factorizable corrections, the non-factorizable corrections 
distort the invariant mass distribution. 



2.1.2 Two-loop order 



At two-loop the separation into factorizable and non-factorizable contributions is somewhat 
more involved. There are four distinct contributions, two of which are non-factorizable. First 
of all there is a two-loop factorizable correction given by the following diagrams 



+ 



+ 



(9) 



Secondly, there are contributions that can be attributed to the interference between one-loop 
factorizable and one-loop non-factorizable contributions. An example of such a contribution in 
the soft approximation is 




+ 




iMc 



k 



2fi 



X 



'^1 '^ii 



h,i2 hh Di + 2/2P1 -l2k2 D2 - 2/2P2 kki -lik'i 



(10) 



The two-loop corrections of this type factorize into a product of one-loop factorizable corrections 
times one-loop non-factorizable corrections. The same factorization holds for hard momenta 
flow in the factorizable loop. Combination of all contributions of this type will lead to 




+ •: 




(11) 
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The first factor, A^i, is tlie same as tfie one-loop factorizable correction. Tlie second factor is 
the same as one-loop non-factorizable correction. 

In addition to the two contributions mentioned above, there are also two classes of non- 
factorizable contributions. The first class consists of all diagrams where two photons are ex- 
changed between the two fermion pairs. Let us list here all these diagrams. There are diagrams 
like 




1 
2 



iXj-^ iXi-^ r\j' 



2^t 



k 



2u 



i,i2 kill kik — /c2^i 



-kol 



2l'2 



(12) 



where the integration is over the virtual photon momenta, li and k-, with the appropriate 
measure, Eq. (|^). We omit factors Di^2/{Di^2 ± '^h'Pi,2 ± 2/2^1,2) in the expression above. 
There are three more contributions of this form corresponding to interactions: {kik'2){kik'2) ■, 
{k[k2){k[k2), {k[k'2){k[k'2). Another set of diagrams contributing to this class is given by 




+ 




1 

'2 



rv-^ r\j-^ r\j' 



2ti 



2v 



kih k'lk -k2li -k2l2 



+ 



ki k\ k2fj, k 



2v 



k'lh kik -k2li -k2l2 



■ (13) 



There are three additional contributions of this type: {kik2){k[k'2), {kik2){kik'2), {k[k2){k[k2). 
Finally, the last set of diagrams contributing to the first class of non-factorizable corrections is 




+ 



^2fi 



2u 



h,i2 k[li kik -k2li ~k'2l2 



(14) 



There is one more diagram of this type: {kik2){k[k2). 

Combining all these contributions we obtain a gauge invariant expression for this part of 
the two-loop non-factorizable correction 



iMo 



ki k[ 



ki k[ 



ko 



ko 



kr. 



h,h \liki lik'i 



J2ki hk'i 
Di 



Do 



-^2^2 ^^2^2 



Di + 2I1P1 + 2/2P1 D2 - 2I1P2 - 2I2P2 
Using the non-factorizable currents this can be written as 



l^Mo [ {Ji-J2){li)*{Ji-J2){l2), 



(15) 




(16) 
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where the currents are the same as in the one-loop case and are given by Eq. (^. The product 
of two currents is defined by the following current multiplication rules 

= jr{h,D)*Jtil2,D) = J^luD) jni2,D±2hp,), (17) 
where ± = + for i = 1 and ± = — for i = 2. Furthermore 

jrih, D) * J^{h, D) = J^{h, D) * Jl^ih, D) = J^h, D) J^{h, D). (18) 



The second class of non-factorizable corrections consists of diagrams with only one photon 
exchange. These are diagrams of the following form 




The loop integration over the fermion loop and over the non-factorizable currents does not 
factorize. 

As we have seen, two-loop non-factorizable corrections receive contributions from two dif- 
ferent sources: from additional interactions with the decay products and from the corrections 
to the propagation of the soft exchange particles. The two contributions given by Eqs. ([l6|) 
and (|1^) respectively are separately gauge invariant and are distinguished in particular by a 
different effective coupling constant. Eq. ( |T6l) originates solely from additional interactions of 
photons with the decay products. This term contains two non-factorizable currents for each 
decay and is proportional to a^. Eq. (|l^) on the other hand receives contributions from propa- 
gator corrections as well as interactions with the decay products. This term contains only one 
non-factorizable current and is proportional to a'^Nf, where Nj is the number of fermions in 
the loop. 

In summary, we separated the factorizable, ~ afact, and non-factorizable corrections, ~ a^f- 
The non-factorizable corrections are further separated into corrections due to the interaction 
of the photons with the (dipole of the) decay products, ~ ctdipoie, and corrections due to the 
propagation of the soft photons, ~ apmp- Thus, the complete two- loop corrections can be 
written schematically as 

= "fact + "factttnf + "nf! Q^nf = "dipole + (o^dipolcaprop) • (20) 

The parenthesis in the last term indicates that the two effects do not factorize. The split up is 
gauge invariant and each contribution has different physical properties. 



2.1.3 N-loop order 



The procedure described above generalizes to the separation of A^-loop correction. The complete 
correction consists of interferences between factorizable and non-factorizable corrections. 
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factorized with respect to each other. Schematically 

N 

«"' = E«f'Ic>M- (21) 

Non-factorizable corrections can be due to the interaction of soft photons with the dipoles of 
the decay products or due to the correction to the propagation of the soft photons 

i 

«M = I](«dipolc«p7op)- (22) 

A useful illustration is given by the four-loop non-factorizable diagrams shown in Fig. |1|. The 
first diagram, Fig. |l|(a), contains four dipole interactions and no propagation corrections. The 
second diagram, Fig. |I|(b), has three dipole interactions and a first-order propagation correction. 
Finally, the last diagram. Fig. |l|(e), has only one dipole interaction and a three loop propagator 
correction. 




Figure 1: Four- loop diagrams contributing to non-factorizable corrections. The di- 
agrams are proportional to (a) ~ Odipolc (b) ~ «dipole«prop, (c) ~ adipole"prop> 
(d) ~ adipoic«prop and (e) ~ adipoieOprop- 



At A^-loop the non-factorizable correction arising solely from dipole interactions, 
can be written in terms of the non-factorizable currents analogously to Eqs. @ and (| 
are given by 



a 



iMc 



1 

m 



[JiJ2]{k) 



N 



N 

dipole ' 

They 



(23) 



(27r)4P 

The products of non-factorizable currents are defined according to Eqs. (|l^) and (|18]). Exami- 
nation of the integral, Eq. (^), confirms that also at higher- loop level the leading contribution 
to the non-factorizable correction is due to soft photon exchange, Q ~ TM/E. 



2.1.4 Charged unstable particle production 

In this subsection we repeat the procedure of separating different effects for the case of charged 
unstable particle production. The difference to the neutral unstable particle case is that now 
in addition to the decay-decay interferences there will be production-decay interferences. 
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Let us consider a toy model with a massive, charged, scalar, unstable particle field, 0. The 
scalar field decays into a pair of massless fermions, ip and if)' . One of the fermions, ip, is charged 
and thus couples to a U{1) gauge field, A^. The other fermion, ■?/;', is neutral. The Lagrangian 
of the model is given by 



IP 



+ 



+ ip + ip' ]^ 
PX(0V) + (p^Viip'^) + h.c. 



(24) 



The difference to Eq. (^ is that now the field cf) is charged, thus p in the kinetic term of the 
(f) field is the covariant derivative. Again, the V and V terms are responsible for the production 
from a neutral source X and for the decay of the unstable particles respectively. As in the 
neutral case, we assume that these terms can be treated perturbatively. 

In order to see how the factorization picture can be generalized to the case of charged 
unstable particles we consider the process 



X ^ 0+(pi)r (P2) ^ ^+{hW{k[)r{k2W{k'^, 



(25) 



where pi^2 = ^1,2 + ^'1 2 ^^'^ P12 ~ ~ TM. For simplicity we will assume that there are no 
charged initial state particles. 

At one-loop the complete non-factorizable corrections can be written as |^ 

dH 



The currents are given by 



iMr 



(27r)4/2 



Jo' 



Pi 



+ 



P2 



for photon emission from the production stage of the process and 



A* 
Pl_ 

IPl 



Iki 



J2' = -e 



P2 



^2 



-lp2 -ih 



D2 - 2lp2 



(26) 



(27) 



(28) 



for photon emission from the decay stages of the process. The main difference is that there 
are three interference terms now. Two terms account for photons connecting the production 
and the two decay sub-stages (production-decay interferences). The third term corresponds 
to photons connecting the two decay sub-stages (decay-decay interferences). The decay-decay 
interferences are similar to those in the neutral unstable particle case. 

At A^-loop order the full correction will be a sum of products of factorizable and non- 
factorizable corrections similar to Eq. (pT)). Non-factorizable corrections come from two sources 
similarly to Eq. (|22|) . One is the interaction of the soft photons with the decay or production 
dipoles. This contribution can be expressed in terms of the non-factorizable currents and is 
given by a generalization of Eq. (^31) . 

1 N 



dH 



(27r)4/2 



[Jo ■ Ji + Jo ■ J2 + Ji ■ J2KI) 



(29) 
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The second source for A^-th order non-factorizable corrections is photon propagation corrections 
to lower-order non-factorizable corrections. Finally, factorizable corrections will be corrections 
to the two decay sub-stages and to the production sub-stage. 

In summary, in the case of production of charged unstable particles the general picture 
still holds. Indeed, at any order the full correction comes from several physically different 
effects. First of all, factorizable and non-factorizable corrections can be separated. Second, 
non-factorizable correction can be production-decay and decay-decay interferences. Each type 
of interference receives corrections from the coupling of soft photons to production and decay 
dipoles and from corrections to the soft photon propagation. 

At this point let us comment on the possibility of a generalization of this separation proce- 
dure to QCD. We believe that in QCD the physical picture should remain the same. At any 
order the corrections should come from the physically different sources identified above and, 
thus, one should be able to separate them. It seems, however, that it will be complicated to 
generalize the procedure we used so far. This is because the separation given above uses ma- 
nipulations of the non-factorizable currents. In QCD the gauge structure is more complicated 
due to the gluon self coupling. At higher order it will be hard to maintain gauge invariance of 
separate contributions. It would be advantageous, therefore, to find a framework, which makes 
use of the physical distinctions between various effects. Such a framework should be manifestly 
gauge invariant. In Sect. § we come back to this question. 



2.2 Integrating over decay products 



In this subsection we study the mechanism that leads to the suppression of non-factorizable 
corrections at high energies, E ^ M. The mechanism is at work for distributions integrated 
over the momenta of the decay products keeping their invariant masses fixed. High energy 
suppression arises because of the coupling of the soft photons to the decay products. Therefore 
we will concentrate on the non-factorizable corrections coming from the interaction with the 



decay dipoles, as defined in Sect. |2.1| . We will not consider any further factorizable corrections 
nor non-factorizable propagation corrections. They are proportional to the coupling constant 
without any extra high energy suppression. This fact illustrates that the effects we identified 
in Sect. ^A] are indeed physically different. Also, we do not discuss high energy behavior of 
production-decay interferences that are present in the charged case. The suppression mechanism 
works differently in this case and we leave this issue for the future. Thus, the estimates that will 
be presented in this subsection are valid for decay-decay non-factorizable interferences coming 
from the interaction with decay dipoles. In Sect. this was shown to be a well defined gauge 
invariant subset of all corrections. 



2.2.1 Neutral Unstable Particles 



The first process we will consider is a neutral, unstable, scalar particle with momentum p 
decaying into a pair of massless, charged fermions with momenta k and k' respectively. The 
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Lagrangian of the model is given in Eq. ([^). We will estimate the A^-th order virtual non- 
factorizable correction due to exchange of iV-photons with momenta li, . . . ,1^. As a first step 
we will integrate the decay part of the Born matrix element squared multiplied by non- 
factorizable currents over k and k', keeping p fixed. According to Eq. (p3D the relevant integral 
is 

r-- sm s(n s\p-t-k'i X Pl^x (A-l-)- . . . {^-^y. (30) 

This tensor corresponds to the decay part of one of the particles. is the Born matrix element 
squared. The indices ai . . .a^ will be contracted with the indices of the non-factorizable current 
from the other unstable particle. 

In the center-of-mass frame at high energy the components of momenta in the formula above 
can be estimated as follows: 

p k k' ^ E, /i ~ TM/E. (31) 

The latter estimate follows from the fact that only soft photons are relevant for non-factorizable 
corrections. Also {kk) ~ {k'li) ~ {pU) ~ TM. Since ultra-relativistic particle decay forward 
k^ and k'^ are almost collinear to and to each other. This will lead to a suppression of the 
non-factorizable currents in Eq. (^). In order to make this suppression explicit, it is convenient 
to change variables 

k^' = ^{p^' + A^"), k'^ = ^{p^' - A"), (32) 
where A ~ M. The phase-space integration over the final state particle momenta is given by 

J d^k d^k' 6{k^) 6{k'^) 6\p-k 
which allows to rewrite Eq. (0) as 

jai---ajv ^ I |-p|2 ( ^ 

J A \kli 

In terms of the new variables, the non-factorizable currents depend on p and A. They can 
be expanded in M/E 

There will be only odd powers of A in the expansion because the exact non-factorizable current 
is anti-symmetric under k ^ k' . 

We are now ready to estimate the size of A^-non-factorizable currents averaged over the 
momenta of decay products, /"i- -"jv^ terms of E, M, and F. We want to compare this to 
the leading contribution of the Born matrix element squared, corresponding to A^ = 0. This is 



-k') = J d^A 5{M^ + A2) 5{2pA) = (33) 



k' \°i / k k' 
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because our goal at the end is to estimate non-factorizable corrections with respect to the Born 
cross section. 



The estimate depends on the parity of A^. If is odd the product of the currents contains 
an odd number of A^^. Upon integration these terms will vanish. For even Eqs. ( ^If ) and 
lead to the following estimate: 

j"i...«iv ^ / °' ^ (36) 
1 Born ■ r~^, for A^ even. 

In Eq. (|36|) we denote the leading contribution to the Born matrix element squared by "Born" . 
Of course, the equation above does not account for the Lorentz structure. However, the infor- 
mation it provides is sufficient to estimate the non-factorizable correction. 

For the purpose of illustration let us make this estimates more precise by taking into account 
the Lorentz structure. As an example, we consider the simplest case N = 2: 



\V\ 



p 



{.Ph){pl2 



(37) 



The first factor is the leading contribution to the Born term while the second factor is the 
leading contribution from the non-factorizable currents. Power counting in the second factor 
gives ~ in agreement with Eq. (|36D. 

The estimate given in Eq. (BO) can be summarized in the following way: using the scaling 



rules given in Eq. (|3T|) the integrals over the decay product momenta, Eqs. (^) and (|3^), can 



naively be estimated to be of the order 0{{E/M ■ l/F)^) relative to the Born term. At high 
energy, however, k is collinear to k', which introduces an extra suppression of 0{M^ / E^) for 
even A^. This results in the final estimate quoted above: 0{T~^) relative to the Born term. 

Finally we are in a position to estimate the non-factorizable correction itself. At A^-th 
order it is given by Eq. (^) and involves A^-integrations over the momenta of virtual particles, 
A^ propagators of soft exchange particles and interference of 2N non-factorizable currents, A^ 
coming from each unstable particle. We assume that one integrates over the momenta of decay 
products. This allows us to use the estimates made in the previous subsection. Denoting the 
non-factorizable correction relative to the Born term by 5nf, we find for even A^ 

The first factor comes from the loop integration measure and A^ virtual particle propagators 
whereas the second and the third factors are due to the 2N non-factorizable currents. The 
non-factorizable correction for odd A^ 

5„f = 0, A^ odd, (39) 

because of the symmetry mentioned above. 
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2.2.2 Charged Unstable Particles 



The estimates made in the previous subsection can be generalized to the case where the unstable 
particle is charged. We will consider a process where the unstable particle, with momentum 
p^, decays into a charged and neutral fermion with momentum and fc'^ respectively. The 
Lagrangian of the model is given in Eq. (p^ . The Born matrix element squared, multiplied by 
non-factorizable currents and integrated over the decay product momenta now reads 

= / - ...(^- -^r. (40) 

J A Kph kliJ \pIn kl^J 

As in the neutral case we change variables and define A'^ through the relation /c = (p + A)/2 
with A ~ M. Expressing the non-factorizable currents in the new variables and expanding in 
M/E results in 



"i^'4^-X (41) 



r V E 

The main difference compared to the case of neutral, unstable particles is that in the charged 
case there are odd powers of A in the expansion. This is because the exact non-factorizable 
current is not anti-symmetric under k ^ k' anymore. As a consequence non-factorizable 
correction will not vanish any longer for odd A^. Instead they will get an extra M/ E suppression. 
The reason for this additional suppression is that one of the non-factorizable currents will have 
to be expanded to second order in A. For even the estimates are the same as in the neutral 
case. Repeating the estimates that lead to Eq. (|36|) one obtains 



Born ■ r-^ ■ M/E, A^ odd. 
Born ■ r^^, A^ even 



N .... (42) 



where "Born" again denotes the leading contribution to the Born matrix element squared. 

As in the neutral case it is possible to get the Lorentz structure. Let us consider, for 
example, the case A^ = 1: 

,2 r 



\v\ 



.2 r 12 



r 



(43) 



{piy 

The first factor is the leading contribution to the Born term whereas the second factor is the 
leading contribution from the non-factorizable currents. Power counting in the second factor 
gives ~ 1/r ■ M/E. 



Using Eq. (|4^) we can now obtain the estimate of the decay-decay non-factorizable correction 
in the charged case. We obtain 

5nf ~ y-^j , A^even, (44) 



5nf ~ (^j , ^odd, (45) 
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These estimates can be summarized as follows: the non-factorizable correction with 2N decay 
currents scales as E'"^^ if is even. Because of symmetry reasons there is an extra suppression 
oi 1/E for each of the two phase-space integrations over the angles if is odd. It should 
be stressed again that this estimate is incomplete. The reason is that we did not consider 
production-decay interferences even though they are present in the charged case. 



3 Effective Field Theory description 

The estimates of non-factorizable corrections presented in the previous subsection seem to be 
quite technical. They show different energy scaling behavior depending on the process, particles 
participating in the process, type of interaction through which they interact and many other 
factors. Moreover, the very definition of what we mean by non-factorizable corrections rests on 



the separation of factorizable and non-factorizable corrections discussed in Sect. |2.1| . There we 
discussed just one specific example (out of many) and showed how such a separation can be 
achieved in QED. 

Despite all the different ways the non-factorizable corrections manifest themselves, all the 
estimates rely on one simple observation. In all processes involving the production of unstable 
particles close to resonance there are always two scales present: a hard scale, Ai ~ E, and a 
soft scale, A2 ~ F ■ M/ E. If E'^ ^ MT there is a hierarchy between these two scales, Ai ^ A2. 
At the hard scale, Ai, the process can be described in terms of interacting stable particles. 
The fact that the particles are unstable is not very important and can be neglected. At the 
soft scale, A2, the instability of the particles becomes important and cannot be neglected any 
longer. If one's objective is to study effects due to the instability then the problem can be 
simplified by separating the two contributions and studying them separately. 

This separation of scales is most naturally reflected in an effective field theory framework. In 
this framework one starts with identifying all the relevant degrees of freedom in the underlying 
theory. They are determined by the scales present in the process. These scales can come from 
two sources. Either they come from the underlying theory itself (masses of heavy particles, 
etc.) or they are introduced through kinematical constraints (CMS energy, etc.). Without 
kinematical scales the mode contribution is determined by the typical excitation energy. How- 
ever, in the presence of kinematically induced scales the excitation energy alone is not sufficient 
to determine the mode contribution. Additional quantum numbers are needed to describe the 
mode contributions uniquely.^ Once the relevant modes have been identified we can integrate 
out the modes that do not show up as external particles. 

Let us recall how the procedure of constructing an effective theory works in the simpler 
case without kinematical constraints. The first step is to integrate out hard momenta. This 



^In the case of heavy quark pair production close to threshold, for example, the kinematical constraints are 
such that particles are produced with a small velocity, v (see |23|). There are three modes distinguished 
by energy with four-momenta: m,^ to), mv,^ mv), rmr,^ mv'^). This is not the complete list, 
however. One needs an additional quantum number — three-momentum — to identify an additional mode with 
four- momentum (~ mv, ~ mv^). 
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results in a theory where only soft degrees of freedom are dynamical. The interactions of the 
soft degrees of freedom are described by an effective Lagrangian 

>^efr = ET^ C»„(soft fields). (46) 

n 

Each operator, On, is multiplied by some power of a scale, A, and a dimensionless Wilson 
coefficient, c„, which is a series in the coupling constant of the original theory c„ = X)c„,iQ;*. 
The contributions of the hard modes are included in the Wilson coefficients. These coefficients 
are obtained by matching the effective theory to the original theory. Since there is a hierarchy 
of scales the operators On are local and A is a typical excitation energy of the hard modes. 

Let us now investigate how this picture applies to processes with unstable particles pro- 
duced close to resonance. This induces a kinematical constraint and consequently the distinc- 
tion between the various modes is more subtle. First of all, modes arc distinguished by their 
momentum, which can either be hard, P ~ Ai, or soft, p ~ A2. By the statement P ~ Ai we 
mean that at least one component of the momentum is of the order Ai and none is much larger. 
However, the fact that there is a kinematical constraint means that an additional quantum 
number is necessary to distinguish different modes. This quantum number is off-shellness or 
virtuality of the particles, D, defined as 

D = P'^- m\ (47) 

where P is the total momentum and m the mass of the particle. Indeed, if the particle is close 
to resonance, P^ ~ m^, it can happen that a particle has large momentum, P ~ Ai, but its 
virtuality is still small, or even zero. In this configuration observables are sensitive to small 
(soft) changes of momentum, P — > P ±p, p ~ A2. The contribution of such modes depends on 
two momenta: the hard on-shell momentum, P, and soft momentum, p. The total momentum 
of such a mode is P ~ Ai, and the virtuality is D ~ A1A2. 

This makes it clear that in order to account correctly for all the relevant modes we will have 
to distinguish between three types of momenta: 

• off-shell hard momenta, P ~ Ai, P^ — ~ A^. 

• on-shell hard momenta, P^ ~ Ai, P^ = 

• soft momenta, p ^ A2. 

In general, the virtuahty can take four values: 0, A2, A1A2, and A^. Consequently, there are 
the following types of modes in the problem: 

• Hard modes: $(P), where P ~ Ai, D ~ A^. 

• Resonant modes $(Pm,p), where P^ ~ Ai, p ~ A2, D ~ A1A2. 

• Soft modes: $(p), where p ~ A2, D ~ A^. 
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• External modes: $(Pm), where Pm ~ Ai, D = 0. 
All fields together with their respective energy and virtuality scales are listed in Fig. ^ 



Ai- 
A2- 




$(P„,p),<l>(P),<l>(P„ 
$(P„,p),<l>(p) 



A1A2 

^*-2 





2 



<|.(P) 

<l>(p) 
$(Pm) 



Figure 2: Relevant energy-momentum and virtuality scales with a list of corresponding 
fields contributing at each scale. 



Let us now consider what happens if we start to integrate out fields in the case where 
kinematical scales are present. It is possible that fields that are integrated out and fields that 
are left in the Lagrangian - being distinguished by the additional quantum number - have the 
same typical energy. In this case the operators On can be nonlocal. The scale A is determined 
by the quantum number, which distinguishes the modes. 

Applying this strategy to processes with unstable particles one first integrates out the hard 
modes, $(P). The resulting effective theory describes instability effects, such as propagation 
of unstable particles close to resonance, $(Pm,p), their decay into on-shell particles, <l>(Pm), 
and interactions between them due to soft exchanges, The effective Lagrangian contains 

nonlocal operators. The main point is that in this way we achieve a separation of factorizable 
and non-factorizable effects. As we will see, the factorizable corrections correspond to the 
hard effects encoded in the Wilson coefficients whereas all non-factorizable corrections are 
described by interactions of the still dynamical soft modes. Since this separation is based on a 
physically relevant hierarchy of scales it is possible to generalize it to more complicated cases. 
The usual way of separating factorizable and non-factorizable corrections is based on ad hoc 
manipulations of various contributions. These manipulations become increasingly tedious for 
more complicated cases. 

As a next step one can integrate out the resonant modes, <l>(Pm,p). This results in a new 
effective theory that contains nonlocal interactions between external particles, $(Pm), coupled 
in the production and decay points. It also contains interactions between production and decay 
points due to soft exchanges, (non-factorizable corrections). Whereas integrating out hard 
modes provides a separation between factorizable and non-factorizable corrections, integrating 
out resonant modes provides a separation of different types of non-factorizable corrections 
(production-decay, decay-decay). 

The last step is to integrate out the soft modes, $(p). This will lead to the 5- matrix defined 
in terms of external particles, $(Pm). 
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3.1 Toy model 



Let us now carry out this programme in some more detail for the toy model specified in Eq. (|l]). 
The list of relevant modes is the following: 

• Hard modes: 0(P), ipiP), A^{P). 

• Resonant modes: (j){PM,p), "ipiPo^p)- 

• Soft modes: ip{p), A^{p). 

• External modes: ipiPo)- 

It should be understood that the ^^(-P) field is subject to gauge transformations at the hard 
scale whereas v4^(p) is subject to gauge transformations at the soft scale. Thus, for example, 
the covariant derivative acting on ip{P) contains A^{P) and the covariant derivative acting on 
ip{P(),p) contains A^{p). 

The leading diagram for the process Eq. (^) is shown in Fig. ^(a) where the dashed line 
denotes a resonant field 0(Pm,p)- Neglecting an overall coupling constant, this diagram is of 
the order (AiA2)"^. 

3.2 Integrating out hard modes 

The Lagrangian given in Eq. ([l|) describes the full theory where all modes are dynamical. The 
first step towards the effective theory we are aiming at is to integrate out the hard modes 
4>{P), i^iP), and Afj,{P). They all have momenta P ~ Ai and virtuality D ~ A^. The resulting 
effective Lagrangian contains fields (j){PM,p) (Pm ~ Ai, P|^ = M^), iI){Pq,p) (Pq ~ Ai, Pq = 0), 
ip{p) and A^{p), with p ~ A2, and the external field ip{Pt))- The Lagrangian of the effective 
theory is determined through matching. Thus, we require that up to a certain order in the 
coupling constants and up to a certain order in A2/A1 the effective theory coincides with the 
underlying theory. At this point we restrict ourselves to contributions that are not suppressed 
by powers of A2/A1. 

Let us start by integrating out 0(P). At leading order in a, we have to consider diagrams 
as shown in Fig. |^. The diagram shown in Fig. ^(a) cannot contribute to the process Eq. (H) 
if the 0-field is hard, since this would violate the kinematic constraint Eq. (^. On the other 
hand, the background process. Fig. 0(b), does contribute to this process, even if the 0-fields are 
hard. Thus, integrating out 0(P) results in an operator X{ipip){il}%p) in the effective Lagrangian. 
However, for each hard 0- field we get a suppression of A2/A1. Thus, in leading order in A2/A1 
there are no new operators introduced by integrating out 0(P). 

Consider now a-corrections in the leading A2/A1 approximation. As shown in Fig. §(a), 
there are diagrams involving hard 0-fields that are not suppressed by A2/A1. However, these 
diagrams do not result in new operators in the effective Lagrangian. Rather, their effect is 
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Figure 3: Resonant (a) and background (b) Born contributions to pair production in 
the underlying theory. Resonant kinematics means that Di ^ D2 ~ A1A2. 



encoded in the (modified) coefficients of already existing operators. In fact, Fig. ^(a) results in 
a modification of VXcpcj). Thus, VXc^xp — > cj^VXcpcj) where c-p = l + 0{a) is a Wilson coefficient. 
Of course, there are also diagrams resulting in a corrections to operators that are suppressed by 
A2/A1. An example of a diagram that gives rise to an a correction to the operator X{il)il))[il)il)) 
is shown in Fig. ^(b). 




Figure 4: Diagrams involving (j){P) that result in (a) a corrections to VXcjxj) and (b) 
a corrections to the operator X{il)ip){'4)ip). In diagram (a) all loops are hard but the 
decaying fields are resonant. In diagram (b) all internal fields are hard. 



As a next step we turn to the U{1) field. In order to integrate out A{P) we have to consider 
diagrams as shown in Fig. ^ For a hard photon, the diagram shown in Fig. ^(a) is suppressed 
by powers of A2/A1. This diagram only gives a leading in A2/A1 contribution for a soft photon. 
The diagram shown in Fig. ^(b), however, is only suppressed by a coupling constant a, but not 
by powers of A2/A1. Thus, it leads to an a correction to the contribution of the (pipVip operator. 
We take this into account by multiplying this operator by a Wilson coefficient cx? = 1 + 0{a). 
Since we are dealing with a neutral field 0, the diagram shown in Fig. ^(c) does not contribute. 

The fields ip{P) enter through self energy insertions of (j){PM,p), thus they will modify 
the operators bilinear in 0. The corresponding operator leading in A2/A1 is 0-P|/0 = 0M^0. 
However, at leading order in a this operator does not contribute. The next-to-leading order 
corrections in a result in an operator (piMTip. This can be seen as follows: the relative order of 
magnitude of a hard self energy insertion is given by aA^A]^^(AiA2)^^ = aAi/A2. The factor 
a ^ comes from the coupling constant, A.\ from the phase-space integration, A^^ from the 
hard fermion propagators and (AiA2)~^ from the additional (f){PM,p) propagator. From the 
estimate above it is evident that the self energy insertion is not suppressed. Therefore, these 
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Figure 5: One-loop radiative corrections in effective theory, (a) would be counted 
as non-factorizable and (b) as factorizable correction in the underlying theory. If the 
unstable particles are charged there are diagrams in the underlying theory that contain 
both factorizable and non-factorizable corrections. An example is given in (c). 



contributions have to be taken into account to all orders. This leads to the usual self energy 
resummation. 

In summary, at leading order in A2/A1 the effective Lagrangian takes the form 
£eff = ^</'(2PmP + ^M^)</) + ^(fg + ji)^ - ^F^,F>^^ + cp</.XP(Pm)</' + cvH^^PkM- (48) 

We should stress that the terms ip j^ijj and F^^F^^'^ in Eq. ( ^Hf ) are not equivalent to the ones in 
Eq. (|l]). In the underlying theory the fields A{P) and ip{P) are dynamical whereas in Eq. (^81) 
they have been integrated out. 

The effective Lagrangian given in Eq. (^8|) is in fact completely equivalent to the double 
pole approximation (DPA) [Q. At Born level this can be seen as follows: The production 
and decay operators, (f)V{PM)<P and cI)iPV{Pm)4', depend on the on-shell momentum Pm- The 
kinetic operator of the unstable particle depends also on the soft momentum p. Consequently, 
the Born process shown in Fig. |^(a) is described by the Lagrangian Eq. ( ^81) as a superposition 
of an on-shell production subprocess, unstable particle propagation subprocess, and on-shell 
decay subprocess. This picture is the same as that in DPA. 

The equivalence of the effective theory given by the Lagrangian Eq. (|48|) and DPA can be 
extended to a-corrections in the leading A2/A1 approximation. In DPA the important fact is 
that a-corrections are separated into a sum of factorizable and non-factorizable corrections. 
The diagram shown in Fig. ^(a) is a non-factorizable correction. Only soft gauge bosons give 
contributions to non-factorizable corrections in DPA. On the other hand the diagram Fig. ^(b) 
is a factorizable correction in DPA. 

Let us see how this happens in our effective theory framework. There are two sources 
of a-corrections. First, the Wilson coefficients contain a-corrections. As we have seen above, 
these correspond to corrections obtained by integrating out hard modes. Second, there are loop 
corrections in the effective theory due to the exchange of soft gauge bosons. In the effective 
theory at leading A2/A1 order gauge bosons couple to resonant particles in the soft photon 
approximation. Thus, the evaluation of a diagram of the type Fig. ^(a) with a soft photon 
is the same in the effective theory and the DPA of the underlying theory. Furthermore, loop 
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integrals in the effective theory have to be carried out in dimensional regularization. A loop 
integral corresponding to the diagram of Fig. ^(b) with a soft photon does not have a scale 
and, therefore, is zero in dimensional regularization. 

All of the above means that the complete factorizable corrections correspond to the contri- 
butions due to the hard modes, calculated in dimensional regularization. These contributions 
are included in the Wilson coefficients of the effective theory. The non-factorizable corrections 
are described by the dynamical degrees of freedom of the effective theory. In this way we achieve 
a separation of factorizable and non-factorizable effects based on the hierarchy of scales in the 
problem. It is clear that the effective theory description provides a separation between factor- 
izable and non-factorizable corrections also at higher orders in a. This is the main advantage 
of the effective theory approach compared to DPA. 

In more complicated cases, for example when the unstable particles are charged, the sep- 
aration between factorizable and non-factorizable corrections cannot be done on the basis of 
diagrams. There are diagrams, like the one shown in Fig. ^(c), which contain both factorizable 
and non-factorizable parts. The separation is still correctly reproduced within the effective 
theory framework. In Appendix ^ we show how the separation works for this example to all 
orders in A2/A1. In order to do that we use an expansion of loop integrals in dimensional 
regularization that is equivalent to effective theory calculations. 

The effective Lagrangian Eq. (^) contains only local operators. On general grounds, how- 
ever, one can expect nonlocal operators to be present. This is because the momentum of the 
fields integrated out, P ~ Ai, is of the same order as the momentum of the fields still present, 
Pm ~ Ai. In order to understand this mismatch let us recall that in our example hard modes 
do not introduce new operators at leading order in A2/A1. Thus, the resonant diagrams are 
described in terms of the same operators in the original and effective theory. These operators 
are local. As mentioned before, hard modes will introduce new operators that are suppressed 
by A2/A1. Consequently, nonlocal operators will appear in the effective Lagrangian only at 
next to leading order in A2/A1. 

This brings us to the question of A2/A1 suppressed operators. As we have shown, the 
effective Lagrangian Eq. (PS]) contains leading operators in A2/A1 and is completely equivalent 
to DPA. The effective field theory approach should enable us - at least in principle - to construct 
A2/A1 suppressed operators order by order in A2/A1. This would open the possibility to go 
beyond DPA and calculate A2/A1 corrections in a regular (in particular, gauge invariant) way. 
It is not the aim of this paper to go beyond the leading in A2/A1 order. However, before 
proceeding, we would like to discuss this point briefly. In the construction of the effective 
Lagrangian we assumed the presence of two widely separated scales. One of the scales, A2, 
was induced by the kinematical constraints that the unstable particles are close to resonance. 
If the unstable particles are far from resonance then there is no extra scale, and the whole 
effective theory approach is not applicable. This means that for different parts of the phase 
space one needs to perform different calculations and switch between the two as one goes from 
one phase-space region to the other. There is a question of accuracy and applicability of the 
calculation in the transition region. This question always arises when effective theories are used 
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for exclusive distributions. If one is interested in inclusive distributions, where one integrates 
over the complete phase space, then the procedure should be the following. The phase space 
should be divided into regions where one of the two calculations is applicable. The calculation 
then is to be performed in the respective regions, assuming extreme scale hierarchy everywhere 
in that region. Then the sum will not depend on the cutoff defining the regions at the order 
at which calculation is performed. It is not possible to perform one calculation, which would 
have the same accuracy in all the different regions of the phase space with a different scale 
structure. In what follows we will be interested only in contributions leading in A2/A1, which 
corresponds to DPA. The possibility of using the effective field theory framework to go beyond 
DPA is interesting and very relevant, but it is still an open problem. 



3.3 Integrating out resonant modes 



Starting from the effective Lagrangian Eq. (^8[) , let us now integrate out resonant modes, 
(f){PM,p), V'(-Po)P)- These modes all have momenta Pm ~ Pq ~ ^i? P ~ ^2, D ~ A1A2. This will 
result in another effective Lagrangian that contains soft modes, ipip) with momenta 

p ~ A2, D ~ A2, and external modes, ipiPo), with momenta Pq ~ Ai, D = 0. The importance 
of operators in the effective Lagrangian is determined by powers of A2/A1A2 = A2/A1. The 
operators in the new effective Lagrangian are nonlocal. The leading operator corresponding to 
pair production of unstable particles is a ?/^^-operator and has the form 

A{x,y,z) ■ ■ XV{y) ■ {^pVij){z). (49) 

This term describes the resonant Born diagram in the leading approximation. The factor V{y) 
corresponds to the production of unstable particles at the point y and {tpVilj){x) corresponds 
to their decay at the point x. As long as the unstable particles are neutral, each of these factors 
is gauge invariant separately. The coefficient A(x, y, z) contains derivatives with respect to x, y 
and z, which we will denote by P^, Py, and P^ respectively. It describes the propagation of 
unstable particles from the production to the decay point and has dimension [—4]. Only the 
virtuality of modes that have been integrated out can appear in the denominator. Thus, A can 
be estimated as 

A{x,y,z) ~ j-^l l + o(^]]. (50) 



(AiA2)n VAi 



As in section |3.2| we will neglect A2/A1 suppressed terms. Of course in this case we know the 
precise form of A. From Eq. ( ^Sf ) we deduce that the propagator of a 0(P/v/,p)-field is given by 
{2PmP + iMT)-\ This entails 

A(a;,2/,^) = 2p^^^ + ,Mr ■ 2P,p, + tMr- ^^^^ 

When the operator P^ acts on a fermion field, ip{P){x), it gives its momentum P, when it acts 
on pair of fermion fields, {il!{Pi)ip{P2)){x), it gives the sum of the momenta. Pi + P2. The total 
momentum of the fermion pair is equal to the momentum of the unstable particle. Therefore 
it has hard and soft components, (Pi + P2) = P + p. Because we integrated out the unstable 
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particle field the effective Lagrangian does not depend on the unstable particle field itself any 
longer. Nevertheless, it can depend of its momentum. 



Of course, in more complicated cases the coefficients in front of the nonlocal operators can 
be quite complicated. An explicit matching calculation is required in order to determine them. 
In Sect. ^.3.1| and Sect. |3.3.2| we write down the nonlocal operators explicitly and show how the 
nonlocal Wilson coefficients can be determined by the matching. Our main goal, however, is 
to reproduce the estimates of non-factorizable corrections presented in Section ^. In order to 
do that we will not need to know the explicit form of the nonlocal coefficients. The estimates 
of various quantities can be determined by the scaling properties alone. In Sect. |3.3.3| we show 
how the scaling properties lead to such estimates. 



3.3.1 Explicit operators and matching in QED 

We start by writing down an effective Lagrangian describing pair production of unstable parti- 
cles. The Lagrangian will be written in terms of (nonlocal) gauge invariant operators dependent 
on soft fields, ipi^p) and v4^(p), and external fields il^{Po)- The soft fields transform locally un- 
der gauge transformations whereas the external fields transform globally. We will perform the 
matching by comparing matrix elements squared calculated in the effective theory and in the 
underlying theory.^ 

The effective Lagrangian consists of kinetic terms and terms describing interactions. It has 
the following structure 

+ A{x,y,z) ■ {^V,^,tlj){x) ■ XVUy) ■ (^Ant^)W, (52) 

where the first three terms are kinetic terms. Note that p contains the covariant derivative 
whereas p contains the normal derivative. Therefore, the soft fermions ipip) do interact with 
the soft gauge field A^{p) but there is no interaction of the external field, ip{Po), with A^{p). 
The last term in Eq. ( ^2] ) describes doubly resonant pro duct ion- decay processes. The field ip 
appearing there is an external field iP{Pq). Of course, one can easily generalize this Lagrangian 
to describe production of more than two unstable particles. The operators "Dint and Vint con- 
tain operators with different powers of the gauge field A^, describing gauge boson radiation 
accompanying the production/decay process. Schematically, they have the following structure 

Vint = v(^l + eF^, + e\F^,y + ..}j (53) 
Ant = v(^l + eF^, + e^{F^,Y + ..}j. (54) 

The first terms correspond to production/decay without radiation. The terms linear in 
correspond to production/decay accompanied by radiation of one gauge boson, and so on. If 

^Matching matrix elements would be incorrect because this does not fully reflect the fact that ^(Pq) are 
external fields by construction. In particular, the color traces have to be taken, because ip{Po) is subject to 
global gauge transformation only. 
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the source X and the unstable particles is neutral then there there can be no radiation from 
the production stage of the process 

Vint = V. (55) 

Let us start with the QED case. We will write down explicitly operators corresponding to 
different terms in ('?/;Dint'i/')(x). The operator corresponding to one photon radiation is 



j dxdydz -X,y-X\z- X) (i){x)Vi){y)) F^^'iz). (56) 



This operator is gauge invariant because ip is an external field subject to a global gauge trans- 
formation and is an abelian gauge field. This operator contributes only to one photon 
radiation. It is convenient to use a Fourier transformed function for the Wilson coefficient 

E(^j(x-X,?/-X|z-X) = / dh e'''i^-x)^is2{y-x)^is^{z-x) s«(si, sslsg). (57) 



The matching calculation in the effective theory and in QED is 
/, a 

k k' 




\ fvt ru L / 



Here we recognize the non-factorizable current structure that we encountered in the QED 
calculations of Sect. |^. The Wilson coefficient can be read off from Eq. (|58D and is given by 



This result generalizes to A^-photon radiation in a straightforward way. The operator respon- 
sible for radiation of photons is 



(V^A„t^)(x) 



J dxdyd^'z S(^i,.,^^,^(x -X,y-X\z,-X,...,ZN-X) (60) 
X (V^(x)r'^(y)) F^i^i(^i) ... F^^'^^(27v). 



This operator is gauge invariant and contributes solely to observable where the number of 
emitted photons is A^. The matching calculation in QED is very similar to the calculations 
presented in Sect. ^ with a product of non-factorizable currents emerging as a result. The 
resulting Wilson coefficient is given by 

(h h'll 1 ^ - ( KiKj_ _ KJ2^\ ( ^mKr, Kn^^n \ |r..^ 

^>.,.,...^.^u^KK,K\L^...Ln} 2^^, yj^^^j^,^^ khk'hl--\klNk'lN klNk'lN)' ^ ^ 

We can summarize the situation as follows: the propagation of soft and external fields is 
described by the kinetic operators in Eq. (^2]). The production and decay of a pair of unstable 
particles accompanied by radiation of an arbitrary number of photons is described by the last 
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term in Eq. (|5^) with the factors ('^/^X'intV') (a;) exphcitly given in Eq. (^). These operators 
reproduce those non-factorizable corrections that are due to interaction with decay dipoles and 
have the non-factorizable current structure encountered in Sect. As mentioned in Sect. ^ there 
are also non-factorizable corrections due to the coupling of soft fermions to soft gauge bosons, 
which induces corrections due to propagation. Since soft fermions as well as soft gauge boson 
fields are still dynamical in the Lagrangian of the effective theory, Eq. (^), these corrections 
are absolutely identical in the effective and underlying theory. 



3.3.2 Explicit operators and matching in QCD 



We now want to generalize the construction of the effective theory to the non-abelian case of 
QCD. The effective Lagrangian has a structure similar to the one encountered in the QED 
case, Eq. (^). However, the explicit form of the operators ('?/^Pint'?/')(a;) has to change now. 
Indeed, the operator responsible for one photon radiation, Eq. (p6|), is not gauge invariant in 



the non-abelian case. This is because the field tensor F^,^ is not gauge invariant if is a non- 
abelian gauge field. In fact it is impossible to construct a gauge invariant operator that contains 
ipipAfj_ if i/j is transforms globally and transforms locally under gauge transformations. This 
is consistent with what is observed in QCD. It merely reflects the fact that the emission of a 
single gluon is forbidden in QCD due to the vanishing of the color trace. The matching equation 
for the emission of a single gluon then is simply 

k 

X*C^^= 0. (62) 




The first non-vanishing process in QCD is the one where two gluons are radiated. Accord- 
ingly, the first non-trivial gauge invariant operator in the effective theory contains two F^^^, 
tensors 



X (iPixWiy)) Tr (^F>'%z,)U{z,, Z2)F''P{z2)U{z2, z^) 
Here we introduced U{zi, Z2), the path-ordered exponential defined as 

U (x, y) = P exp 



(63) 



ig / Ai,{uj)duj^ 



(64) 



in order to ensure the gauge invariance of the operator. The matching calculation in the effective 
theory and in QCD for the two gluon process is similar to the one in QED 



Tr (W) S(1^(A;, k%, h) r.g"" - 1^,9^(129'^ - 1'29''^) + (1 - 2) 
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The solution for the Wilson coefficient can be read off from Eq. (|65| ) 

^(2) n- (J^iX_ KA_\ ( k^K KK \ , . 



This is similar to what we found in QED. However, there is an important difference compared 
to the QED case in that now the operator that is responsible for the emission of two gluons, 
Eq. (^), will contribute also to processes with a larger number of gluons radiated off. To 
illustrate this point, let us consider a process with the emission of three gluons. The diagrams 
of the effective theory contributing to this process are shown in Fig. |^. There will be three 
distinct contributions. In the first one, shown in Fig. H(a), two gluons are emitted due to 





(a) ■ (b) ' (c) 

Figure 6: Diagrams of the effective theory contributing to the three gluon process. 



the {ijj'ip)F^ operator, Eq. (|63D, with one of the gluons subsequently splitting into two due to 
three-gluon interaction coming from the kinetic term of the gauge field, —1/4 F^j.F'^^. The 
second is shown in Fig. |^(b) and is also due to the {ipil))F'^ operator. There are two sources 
in {%l>'il))F'^ that give rise to terms contributing to three gluon radiation. First, the non-abelian 
commutator terms in F^j^ result in terms with three (and four) gauge fields. Second, the 
expansion of the path-ordered exponentials ^7(2:1,^2) creates terms with an arbitrarily large 
number of fields. Thus, the operator (■?/''?/') given in Eq. ( |63|) contributes to all processes 
with two or more gluons emitted. Finally, the third contribution to the three gluon process 
is shown schematically in Fig. ^(c). This contribution is due to a new ipipF^ operator with a 
Wilson coefficient T,^^^ defined as 



= / dxdydh Sg),^^^,^^3,3(x -X,y- X\z, -X,Z2- X, - X) (67) 
X U{x)V^l,{y)) Tr (f^-''-{z^)U{z^, Z2)F>^^'''\z2)U{z2, z^)F^''''{z^)U{z^, z^)\. 



(V^AntV^)(X) 



The sum of the three contributions shown in Fig. |^ should be equivalent to the full QCD 
calculation. This matching condition determines the Wilson coefficient Tj^^\ 

This picture is the non-abelian version of what we found in QED. Non-factorizable correc- 
tions are either due to the interaction of soft gluons with the decay dipoles or due to propagation 
corrections of the soft gluons. In the abelian case the propagation corrections arise solely from 
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interaction with soft fermions. In the non-abehan case there is an additional source of propa- 
gation corrections, namely the self-couplings of gluons. What we achieved is a generalization 
to the non-abelian case of a gauge invariant separation between non-factorizable corrections 
due to interaction with decay dipole and non-factorizable corrections due to propagation. The 
analogy of the estimates of Sect. ^ for the non-factorizable corrections due to the interaction 
with decay dipoles is now simply an estimate for the Wilson coefficients T.^^^ For example, 
the contribution shown in Fig. ^(c) is due to a 0{g^) dipole interaction. On the other hand, 
the contributions shown in Fig. |Ka) and (b) are due to Oi^g"^) dipole interactions plus 0{g) 
propagation corrections. Thus, it is only within the effective theory framework that the true 
meaning of the estimates for non-factorizable corrections in the non-abelian case becomes clear. 

The matching calculation in QCD is far more complicated than in QED. What will be 
important for us here is that the Wilson coefficient can be estimated as 

^SliM2.2M3.^3(^' ^'1^1' ^2, h) ~ 9 %. \ ^6 ■ (68) 



In the following section we concentrate on the estimates of the non-factorizable corrections 
due to interactions with the decay dipoles. It will be sufficient to know the estimates of the 
Wilson coefficients, Eq. (68), as well as the Lorentz structure and the form of the operators. The 
explicit form of the matching coefficients will not be important. Nevertheless, we performed an 
explicit matching calculation for S^^^^ in order to convince the reader that the effective theory 
does exist and can be matched to QCD. We should mention that the set of operators we have 
chosen is sufficient only at leading order in A2/A1. If we were to go beyond this approximation 
additional operators with a more complicated 7-matrix structure would be needed. Details of 
the matching calculation are given in Appendix |B[ 



3.3.3 Scaling properties and estimates 

In this subsection we will generalize the estimates of the Wilson coefficients found above and 
then apply these in order to obtain estimates of the non-factorizable correction. Thereby we 
will confirm the results that have been found in Section |] by analyzing Feynman diagrams. 

Instead of working with the momenta of the fermions it will turn out to be more convenient 
to work with the momenta Pm + P and A, defined as follows: 

Pm+P = Pi + P2. A = k-k', PmA = -pA, (69) 

where Pm is the hard on-shell momentum of the resonant 0-field that decays into the fermions 
with momentum k and k'. p is the corresponding soft momentum, p ~ A2. Finally, A is a hard 
off-shell momentum, A ~ Ai, however PmA ~ A1A2. 

The estimate found in Eq. (pS]) can be generalized to the Wilson coefficient S*^"-* of the 
operator [il)il))F"' 

~ ^"(AiA2)2" ^t^iKi ■ ■ ■ f^fj-sKs ~ ]y2^ ■ ^^^^ 
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Rewriting this in terms of the new momenta we find 



{Pm or A)2"+>"^ 
(AiA2)2"+™ 



(71) 



We would hke to stress again here that this selection rule is valid only up to A2/A1 suppressed 
terms. 

The selection rule, Eq. ([7T|) , can be understood as follows: the amplitude for radiating one 
soft gauge boson is proportional to gPf^J {Pmv) ~ 9/^2, plus A2/A1 suppressed contributions. 
Consequently, the Feynman rules corresponding to such an interaction operator should scale as 
~ This means that the term linear in gA^j_ should scale as gAj K2. As is a soft field, 

this is equivalent to terms linear in gP^v scaling as gP/ K\. These arguments can be extended 
to radiation of ra-gauge bosons. This leads to the observation that {gP)'"'-teim.s have to scale 
as {gP)'"" / as stated in Eq. ([70|) . Taking into account that only A1A2 can appear in the 
denominator we balance the dimension of the contribution of the operators by powers of Ai 
and A2 in the numerator. Powers of Ai correspond to hard momenta P or A whereas powers 
of A2 correspond to soft momenta p. This leads directly to the selection rule Eq. ([7T|) . 

Let us use this selection rule to study the non-factorizable correction due to the exchange 
of one soft field. As mentioned above, Eq. (|62D, this contribution will vanish in QCD and 
our estimate will only be relevant for QED. The Feynman diagram corresponding to this one- 
loop correction is shown in Fig. 0(a). According to Eq. (|7TD, the leading contribution to the 
Feynman rule for the vertex scales as 



(J'm.P.A) 



A2 



More explicitly, the following operators linear in P^^, are possible 

{^V,^t^){PM,p, A) 



A1A2 



+ 



(A/^p- + P^,pn Pi + A2 



A1A2 



A1A2 



+ 



AiA: 



(72) 



itpVij)P,,, (73) 



where in the last term we used the fact that PmA is subleading. The leading contribution to 
the Feynman rule is coming from the first terms in Eq. ([73|). With the help of Eq. (|72D the 
one-loop correction relative to the Born term can be estimated as 



(i^/ e e 
P T2 A2 



. „ . „ „ A2 e e 
^^^' — -- ^AlA^A^^' 



(74) 



where / is the soft momentum of the exchanged gauge boson. This shows that non-factorizable 
corrections to exclusive quantities are 0{a) and are not suppressed by any additional parame- 
ters. 

In order to compare to our previous results, let us recall that the estimates of non-factorizable 
corrections presented in Sect. || are valid for a special class of observables only. The estimates 
are valid if one integrates out part of the phase space corresponding to final the state fermions, 
keeping the momenta and invariant masses of the unstable particles fixed. This corresponds to 
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Figure 7: Feynman diagrams for one-loop (a) and two-loop (b) non-factorizable cor- 
rections in an effective theory where resonant modes have been integrated out. 



the construction of operators in the effective field theory using Pj^j and p^, but not A'^. Indeed, 
suppose that the contribution of an operator depends on A. After performing a phase-space 
integral over the fermion decay angles this contribution will be expressed in terms of Pm and 
p only. Thus, in order to reproduce the estimates of Sect. Q we will have to construct a new 
effective Lagrangian. This Lagrangian is the same as the previous one, but without terms con- 
taining A. This Lagrangian will be used to calculate observables, in which the integration over 
fermion decay angles is already performed. Consulting Eq. (|73|) , the leading operator linear in 



P^lu is 



(V5Pint^)(PM,p) 



i^pVtP) F^,. (75) 



^2) 

Therefore, the Feynman rule corresponding to the one-gauge-boson-radiation vertex in this new 
Lagrangian scales as 



{Pm,p) ^Ve 



1 P2 



M 



A2 A^ 

The estimate of the non-factorizable correction in this case is 

' Ai As Aj A2 Af E^- ^ ' 

This is an estimate similar to what we had in Sect. |^. It shows that at high energy, E ^ M, 
non-factorizable corrections are suppressed. 

Of course, it can be that the underlying theory has more symmetries. In that case one has 
to take them into account in the construction of the relevant operators. This is the reason 
why the estimate obtained in Eq. ([75]) does not agree with Eqs. ( PBD and (pPj). Indeed, if 
the unstable particle is neutral and it decays into a fermion-antifermion pair with the same 
flavour there is an antisymmetry under the exchange of final state particles. This leads to 
the statement that the one-loop non-factorizable correction is zero in this case. In fact, all 
odd-order-loop non-factorizable corrections are zero because of this symmetry. In the effective 
theory language this means that operators linear in F, like the one given in Eq. ( |75D are forbiden. 
Alternatively, performing the matching one will find that the corresponding Wilson coefficients 
vanish. However, if there is no additional symmetry the estimate given in Eq. (ITS) holds. 
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Let us extend this analysis to two-loop corrections. The explicit form of the terms quadratic 



in Ff,^ is 



= e' jj^, ( 9'"'pI:,p^,p'm + g'^'g^'Ptf ) m f,.f^,. (77) 



The diagram corresponding to two- loop correction is shown in Fig.|^(b). The Feynman rule for 
the vertex with the emission of two gauge bosons scales like 




p2 p2 , p2 



(P„.p) ~ P^.^^l + ^j. (78) 

At high energy the leading contribution comes from the g^" P^Pj^Pl^ term and the contribu- 
tion of the g'^^g^PPf^ term is suppressed by M'^/E'^. The relative two-loop non-factorizable 
correction can be estimated as 

d^d^ f^ph\ (^Ek\ rMV (^Ek\ (^Ek\ .4 .7q^ 

kl kl Ui Af AaI Af y Ui; \Kl Kl ) \kl Kl ) ^ E^' ^ ' 

This is precisely the estimate for two-loop non-factorizable final-final corrections obtained in 
Sect. I, Eq. (PD. 



It is possible to construct A^th order terms in F^j^ for Pint- It is clear that at high energy the 
leading contribution will come from operators in which as many P^,j as possible are contracted 
with Ff^i^. This is because the remaining P^^ needed for dimensional reasons will have to be 
contracted with each other. Since P^j ~ <C A^, this results in a suppression at high energy. 
If is even {F^y)^ has 2N indices to be contracted. Due to the antisymmetry of the field 
tensor a maximum of A^ of these indices can be contracted with P'^j. The remaining P^^ have 
to be contracted with each other. By looking at the selection rule, Eq. (|7lD , it can be seen that 
different m are possible. All of them give the same contribution if is contracted with P'^j. 
If is contracted with F^^, then some extra would have to be contracted with each other 
giving additional M/ E suppression. This means that the leading correction is the same as that 
coming from the selection rule with m = 



(^Pint^)(PM,p) 



p2N 

^ {iljV^){eFf, N even. (80) 



\2N 



Note that this is in perfect agreement with the explicit result of S^^-* given in Appendix If A^ 
is odd (P^i/)^ has 2N indices to be contracted, out of which a maximum of A^ can be contracted 
with Pj^. The rest of the indices cannot be contracted with each other because their number 
is odd. Thus, at least one p^ is necessary. Again, in the selection rule Eq. ([71|) different values 
of m are possible. However, the leading contribution comes from m = 1. All other values of m 
give the same result or suppressed contributions. The selection rule Eq. ([7T|) gives for odd A^ 



(V'Ant^)(PM,p) 



P 

(AiA2)2^+i 



{^/jV^){eFf, N odd. 
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Thus, the operator V-^^t is of the following form 

(^Pint^)(PM,p) = 

= {^^JV^p) + e 



12) 



+ e 



N 



(AiA2)= 



JDl^ TDV p2 

2 \Mf_Mf_M_ 

(AiA2)^ 



M 



(AiA2)2^ 



T 



cri...o-jv 



Filial- ■■F^ii^o-]^ 



N=even 



+ e 



N 



(p2 \(Af+l)/2p/ii pMJV , 
(AlA2)2^+l 



+ 



N=odd 



where To-i...(jjv ^ symmetric tensor of even rank constructed from (y'o-io-, - 



We are now ready to estimate the A^-loop non-factorizable correction. The estimate will be 
different for odd and even A^. Let us start with even N. The Feynman rule corresponding to 
the radiation of gauge bosons from the decay dipole is 



N 



Correspondingly, the relative A^-loop non-factorizable correction is 



N even. 



^3) 



kl ■ 

A|^^ 



A-2 



(Pm 



^{PM^ 

A^ Af 

2 ^JV/2^ 2 



A^ 



A^ 



N 



^2N 



A^ 

\e) 



Af 



N 



N even. 



(84) 



Here the first factor comes from the A^ loop integrals and A^ gauge boson propagators while 
the second factor comes from two interfering currents to radiate A^ gauge bosons. 

The case of odd A^ is treated analogously. The Feynman rule corresponding to radiation of 
N gauge bosons from decay is 



A^ 



V2 Af+1 

In this case, the relative A^-loop non-factorizable correction is 



N odd. 



^5) 



6 



N 



(£h d^kN f e~ {Pj 



,2^)(A.+l)/2, 



^2 



Va^ 



Af+l 



A^ 



A|^^ 

A2 



-.N 



A^ 



Ar 



E J 



N odd. 



^6) 



Thus, the estimates obtained using the effective theory, Eqs. ( plf ) and (|86|), agree with the 
estimates obtained through the analysis of higher order Feynman diagrams, Eqs. (|^) and (^5]). 
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4 Applications 



In this section we consider some topical applications of the techniques presented in Sect. ^ and 
Sect. 1^ to the pair production of W and Z bosons and top quarks in 77 and e~^e~ collisions.^ 
As above, all our considerations will be performed in the pole-scheme. We shall also confine 
ourselves to the case when no kinematic restrictions are imposed on the configurations of the 
final state particles. In general, as shown in 0, the kinematic cuts applied to the final state 
particles may induce significant changes in the energy behavior of the non-factorizable effects. 
We focus on the distributions in which the integration is performed over the momenta of the 
decay products keeping tracks of the invariant masses of the unstable particles. 

In Sect. and Sect. |^ it was shown that non-factorizable corrections due to interaction 
with decay dipoles are suppressed at high energy. Throughout this section we consider only 
this type of corrections. We do not discuss other corrections (factorizable corrections and non- 
factorizable corrections due to propagation corrections) for which the suppression mechanism 
is not applicable. 

4.1 Neutral unstable particles 

We begin with the basic case when only neutral (both with respect to electromagnetic and 
color charges) particles participate in the production process and, therefore, only interference 
between radiation accompanying the decays occurs. This can be exemplified by the process 

77 — > Z{pi)Z{p2) — > 4 fermions (fci, k[; ^2, K)^ (S"^) 

with pi 2 = /ci,2 + ^12- '^^^ non-factorizable effects appear in the first order in QED and in the 
second order in QCD. In the case of color singlet unstable particles the QCD interconnection is 
additionally suppressed by l/N^ compared to the leading QCD contribution, see e.g. |18| {Nc 
being the number of colors). 

For the purposes of illustration it is sufficient to focus attention on the pure QED phenomena 
accompanying the process 

77 ^ Zip,)Z{p2) ^ e-^{h)e-ik',)fi+ih)fi-{k',). (88) 

We present below an estimate of the energy behavior of the non-factorizable correction of order 
0{a^) induced by the exchanges of N virtual photons with momenta between the 

decay products of the two Z's. 

Similarly to what we did in Sect. |^ in order to find the A^-loop contribution we have to 
evaluate the integrals over the Born decay matrix element squared multiplied by the non- 
factorizable currents corresponding to radiation off the massless charged fermions with momenta 

^Our approach may be generalized to the case of pair production of unstable particles in hadronic reactions 
(for discussion of the gluon bremsstrahlung effects accompanying the tt production at hadron colliders see for 
instance [Q). 
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ki 2 and k[ 2 keeping the momenta of the Z bosons, pi 2, fixed 

= J d% d% 6{k^) 5{kf) 5\p, -h- k',) X A^, 



kih k'AiJ \kilN kU 



In 



19) 



where i = 1,2. 



The new element with respect to the consideration of Sect. 0is the appearance of the decay 
tensor A^j,, which resuhs from the vector nature of the decaying particles. The indices fi, v 
will be contracted with the rest of the Born matrix element squared. The indices a\. . .a^ will 
be contracted with the corresponding indices of the non-factorizable currents from the second 
unstable particle. A^,^ can be decomposed into two contributions with different P-parity 

A,. = Aj-, + (90) 

^lu ~ KK + - {kk')g,,, ~ e,,^,k^k'P. (91) 

The A^^, piece corresponds to the iV'^ + A^) contribution to the matrix element squared while 
A^j, results from the (^^4) interference. We shall call them {V'^ + A^) and {VA) contributions 
respectively. 

Similarly to Sect, ^we introduce the new variables 

K = \{pt + ^^\ k't = \{P^ - ^^)- (92) 



In terms of these new variables the decay tensor components can be rewritten as 



A-=l 
2 



A^^ = -6,,.,AV. (93) 



1 

— ( 

2 

We immediately see from Eq. (|93| ) that the A-dependent part of the {V"^ + A'^) component is 
suppressed by a factor ~ M'^/E'^ as compared to the leading term while the {VA) component 
is suppressed by 0{M/E) factor with respect to the leading (V"^ + A"^) part. 

Let us perform an integration over the final state particle momenta in Eq. ( ^9]) (see Sect. 
and, in particular, Eq. (PO) and Eq. (P^)) and evaluate the high-energy behavior of different 



components of the A^-loop tensor 1^^"'°''^ 



ViA)ja....a, ^ f ^ViA) . (k_ _ f±_ _ 

A \kh k'hJ '"\klN k'lNJ ■ ^ ^ 



The result depends on whether N is odd or even 



^ 1 Born ■ T"^, for even. 



(95) 



Born ■ r-^ ■ M/E, for odd, 
''^'^ 1 0, for even. 



^uy 1 n f^^ AT „,.„v. \-'^) 
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Here "Bom" is a symbolic notation for the leading {V^ + A^) contribution to the Born matrix 
element squared. 

Power-counting rules established in Sect. ^ allow one to estimate the suppression factors 
corresponding to the different parts of the non-factorizable correction to the cross section of 
the process (^Sj). The results are summarized in Table where the notations {V^ + A^) - 
{V^ + A^), [V"^ + A^) - {VA), {VA) - (VA) are used to specify the interference between the 
different contributions into the Born decay matrix element squared of each Z-boson. 



N 


Born 




(V^ + A^) - (VA) 


{VA) - (VA) 


A^-even 










A^-odd 








{M/Ey^+'' 



Table 1: Estimates of the energy dependence of the 0{a^) non-factorizable correc- 
tion, 5nf) integrated over the momenta of the decay products for different interference 
contributions corresponding to the process (p^). 



As follows from Table |l|, the one-loop QED non-factorizable corrections to the cross section 
of the process (|88|) for different final states may appear only as a result of the (VA) - (VA) 
decay interference and in practice are numerically small (see Ref. [|16|)- This term scales as 
with the CMS energy, which is in accord with the earlier observation in Ref. . The two-loop 



QED interconnection effects are caused by the interference of the P-even terms [V"^ + A"^) - 
[V"^ + A^), and their contribution also scales as E~'^. Symbolically, the four leading QED terms 
in relative non-factorizable correction, 5nf, for the process can be presented as 



5, 



nf 



ta a 



M 



4 r 



E J 



EJ 



+ ... 



+ a' 



My 
'eJ 



aH-) +... 



(97) 



where parameter is determined by the product of vector and axial couplings of the cor- 
responding neutral current. Therefore, the actual parameter of the perturbation theory in 
5nf separately for even and odd powers of a is a'^{M/E)*, and the series in Eq. ( P7| ) rapidly 
converges at E ^ M. 

All these observations remain valid for the decay-decay part of the QED non-factorizable 
corrections to the process 



e e 



Z(pi)Z(p2) 4 fermions (ki, k[; /c2, 



(91 



The difference is that now there are also production-decay QED interferences, which were 
absent in the process (pS)). Estimates for production-decay interferences are in general different 
from those for the decay-decay interferences. 
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Let us now turn to the QCD non-factorizable corrections to the processes (^) and (|98D 
with final state quarks. Due to the group structure of QCD these can appear only at > 2, 
see e.g. 0]. Only the {V'^ + A^) - {V"^ + A^) decay-decay gluon interference survives at 



N = 2 with the energy scaling E . As we discussed in Sect. |3.3.2| , in QCD the separation 



between higher-order non-factorizable corrections due to interaction with decay dipoles, which 
is of interest to us here, and non-factorizable propagation corrections is complicated due to 
the non-abelian structure of the gauge group. This separation is correctly understood in terms 
of different operators in the Lagrangian of the effective theory discussed in Sect. p.3.2| . What 
we estimate here are the leading non-factorizable corrections due to operators describing decay 
dipole interactions. Using Table |l| we can see that the four leading terms in S^f''^ can be written 
symbolically as 



M 

'e 



/MY 



31 M 



M 



l + al{--] +... 



4 



(99) 



The Nc dependence of non-factorizable corrections might be of special interest. The exact 
formula is complicated for arbitrary order, N, however, the leading Nc dependence is easy to 
estimate and goes as [Cfois)^ /N"^ ~ {Nca^)^ /N"^. This means that the actual parameter of 
the perturbation theory separately for even and odd powers of as is N^al{M/Ey. The same 
conclusion remains valid for the QCD non-factorizable corrections to the process 

e+e-[77] W+W~ 4 quarks. (100) 

The main qualitative difference between the W^W~ and ZZ interconnection results stems from 
the difference in the vector and axial couplings of the W and Z bosons. Note that since the 
Z mass and other properties are well known, ZZ production provides a unique laboratory for 
studying QCD interconnection phenomena |TE|, pi| . 

Formally, the application of perturbation theory to the description of non-factorizable QCD 
phenomena accompanying ZZ and W~^W~ production at ^ M is restricted by the require- 
ment that the typical space-time separation between the decay vertices does not exceed the 
characteristic strong interaction scale i?ch ~ l//^? where ~ 171^, ~ ^qcd, 

d ~ 1/TE/M < (101) 

see Refs. P, |18|. In other words, remembering that the life-time of an unstable particle in the 
laboratory frame is 

El , , 

tdec ~ j^Y^, (102) 

we can say that this condition implies that the decay happens before the formation of the first 
interconnecting hadrons. 

At the same time, since the relevant momentum scale for the non-factorizable effects is 
k ^ 1/d the requirement Eq. ( |101| ) coincides with the Landau pole condition for as associated 
with the S^f'^ correction. Therefore, even though formally the perturbative results at very 
high energies are not justified our estimates show that in the non-controllable domain the 5^*"^ 
corrections due to interaction with products of decay are decreasing very sharply with increasing 
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energy, if the two outgoing hadronic systems are more and more boosted apart. Moreover, the 
higher is the power of as couphng to the decay dipoles the steeper is the energy fall-off. The 
corrections to propagation, on the other hand, are not suppressed at high energy. At large as 
the description of propagation becomes purely non-perturbative. In is worthwhile to mention 
that the non-perturbative models predict that the hadronic cross-talk between the decaying 



particles dampens with energy comparatively slowly 19 



4.2 Charged unstable particles 

In the case of charged, unstable particles both types of non-factorizable corrections are present: 
decay-decay and production-decay. We shall consider only the decay-decay part. Our analysis 
embraces, in particular, QED effects in the W^W~ production process and QCD effects in the 
production and decay of a ti pair 

ti bW^bW-. (103) 

Let us consider the case when the "charged" (with respect to QED or QCD interactions as 
appropriate) unstable particle with momentum decays into a "charged" massless particle 
with momentum k^^ and a "neutral" particle with momentum k'^. We begin with vector-boson 
pair production accompanied by the leptonic decays. Analogously to the previous subsection 
we introduce the new variable A (see Eq. (^2])) and evaluate the high-energy behavior of the 
A^-loop tensors ^(^)jp-"iv 

Ja \ph klj \plN kl^J ^ ^ 

Here the contributions with different P-parity are determined by the Born decay matrix el- 
ements. Repeating the same procedure as in the derivation of Eqs. ( p5| ) and (^) we arrive 
at 

V „^ r Born ■ r-^ ■ M/E, for odd N, 

t"" ~\Born-r-^, for even AT. ^^^^^ 

Aja,...a^ r Born ■ r-^ ■ M/E, for odd A^, 

^ \ Born ■ ■ MVE2, for even A^. ^ ^ 

Again, "Born" is a symbolic notation for the leading {V"^ + A^) contribution to the Born matrix 
element squared. Recall that the main difference with respect to the case of neutral particles is 
caused by the presence of the odd powers of A in the expansion of the decay currents over A. 

The estimates for the energy scaling of the non-factorizable corrections caused by the final- 
final interferences are summarized in Table ^ 

From Table ^ we can see immediately that the leading part of the decay-decay interference 
contribution to 5^^^ in the case of the W^W~ production scales with energy as E""^. This 



was first established in Ref. [|17|- Symbolically the four leading QED terms for this process can 
be written as 



)nf ~ a 
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4r n/M\4 
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My 
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1+a^f-) 



(107) 
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Born 




{V'^ + A^) - (VA) 


iVA) - (VA) 


N-even 






{M/EY^+^ 


iV-odd 




{M/Ey'^+'' 


{M/E'f^+-' 



Table 2: Estimates of the energy dependence of the 0{a^) non-factorizable correc- 
tions, integrated over the momenta of the decay products for the decay-decay interfer- 
ence contributions associated with the "charged" unstable particle production. 



A similar argumentation remains valid for the QCD non-factorizable corrections correspond- 
ing to the ti production process 

77(e+e") ^ tt ^ bW+bW, (108) 

with a substituted by as- Recall that in the ti case there is no color suppression in if 
the pair is produced in the color singlet state and there is an additional suppression for 

the case of the octet ti state [0, |T5[. Analogously to the case of color-singlet unstable particle 
production discussed above the perturbative results for S^[^^ are fully controllable only if the 
restriction Eq. ( |101| ) is satisfied. 



We stress again that for the charged particles case this estimates are not complete. For the 
evaluation of the magnitude of the total QCD interconnection effects one has to incorporate 
also production-decay interferences and real radiation. 



5 Conclusions 

In the near future, the direct evaluation of two-loop corrections to processes involving un- 
stable particles is not feasible, not to mention even higher order corrections. In view of the 
phenomenological importance of these processes, however, it is essential to improve our under- 
standing beyond the presently available one-loop calculations. In this paper we took a first step 
in this direction by analyzing higher-order radiative corrections to pair production of unstable 
particles with their subsequent decay. 

Having a future linear collider in mind, we are particularly interested in the high-energy 
behavior, E ^ M, of the corrections. We identified various effects and studied their behav- 
ior/suppression at high energy. In particular, we considered distributions that are obtained by 
integrating over the angles of the decay products, but keeping their invariant mass fixed. 

As a first case, we studied neutral, unstable particles decaying into fermions that interact 
through the exchange of abelian gauge bosons. This is the simplest case and can be analyzed 
by studying Feynman diagrams. By extending the well known one-loop analyses 
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we found that to all orders in perturbation theory the complete corrections have the following 
structure: there are factorizable and non-factorizable corrections and interferences between 



them that factorize with respect to each other, Eq. (0). The factorizable corrections are 
associated with a hard scale ~ E, whereas the non-factorizable corrections are associated with 
a soft scale ~ TM/E. Furthermore, there are two sources of non-factorizable corrections, 
namely interactions of the gauge bosons with the decay dipoles and propagator corrections of 
the gauge bosons. 

Since the various corrections listed above have a different physical origin it is to be expected 
that their high-energy behavior is different. Indeed, we found that the corrections due to 
the interaction with the dipoles is — on top of the usual coupling constant suppression — 
additionally suppressed at high energies, Eqs. (|38|) and (p9|). The other corrections identified 
above do not have such an additional suppression. This analysis can be generalized to the case 
of charged, unstable particles and similar results are found, Eqs. (H) and (^51). 



The extension of such an analysis to the non-abelian case is a rather daunting task. The 
reason is that the self-coupling of the gauge bosons complicates the matter enormously. In 
particular, most diagrams will now contribute to several types of corrections, rather than just 
one. Therefore, it was necessary to tackle this problem in a different way. Since the separation 
of the full corrections into various parts relies on the presence of two different scales, Ai ~ E, 
A2 ~ TM/E with Ai ^ A2, an effective theory approach looks very promising. 

In order to pursue this we identified the relevant fields present in the underlying theory: 
hard, resonant, soft and external. By subsequently integrating out the various fields we obtain 
a hierarchy of effective field theories. The first step is to integrate out the hard modes. This 
leaves us with a theory where only resonant, soft and external particles are dynamical. The fac- 
torizable corrections are those that are encoded in the (higher order corrections to the) Wilson 
coefficients, whereas the non-factorizable corrections are the corrections due to the still dynam- 
ical degrees of freedom. What we have gained is that the generalization of this separation to 
higher orders in the coupling constant and to the non-abelian case is straightforward. Through- 
out the paper we restrict ourselves to the leading terms in A2/A1. In this approximation, the 
effective theory is equivalent to the DPA of the underlying theory. 

The second step in constructing an effective field theory is to integrate out the resonant 
modes. The corresponding effective Lagrangian contains again (non-local) gauge-invariant op- 
erators, multiplied by Wilson coefficients. Some of these operators are responsible for the 
propagator corrections and some are responsible for the dipole interactions. Thus, in the same 
way as integrating out hard modes provided us with a separation of factorizable and non- 
factorizable interactions, integrating out resonant modes allows us to separate the two kinds 
of non-factorizable corrections. As a consequence of the gauge invariance of the operators, this 
separation is also gauge invariant. 

In order to obtain estimates of the various corrections similar to those obtained by analyzing 
Feynman diagrams, all we have to do is to get estimates of the Wilson coefficients of the cor- 
responding operators. These estimates can be obtained by general considerations, i.e. without 
explicitly knowing the Wilson coefficient. Thereby we generalize the estimates obtained in the 
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abelian case to the non-abelian case. Obviously, if we wanted to do a precise calculation in 
the effective theory we would need the exact form of the Wilson coefficient. They have to be 
obtained through matching the effective theory to the underlying theory. We have done this 
matching for the operators corresponding to the emission of up to three gauge bosons, thereby 
confirming the estimates obtained by general considerations. 

We applied the estimates to several phenomenologically relevant processes. The actual 
estimate depends on the structure of the Born term. Also, it is important to stress again that 
the estimates obtained here are valid only for the non-factorizable corrections due to dipole 
interactions for distributions integrated over the decay angles. Furthermore, in the effective 
theory approach we considered the case of neutral, unstable particles only. 

In the case of charged, unstable particles there are additional non-factorizable corrections 
of the production-decay type. In the effective theory language this means that the operator 
associated with P, Eq. (|53|), is more complicated than given in Eq. (|55D. Allowing for the emis- 



sion of gauge bosons from the unstable particle production process results in more complicated 
operators in the effective Lagrangian. In order to obtain estimates in this case one would have 
to construct these operators and obtain estimates of the corresponding Wilson coefficients. 

Another possible future development is the inclusion of real gauge-boson radiation. As 
mentioned in the introduction, real gauge-boson radiation contributes to the non-factorizable 
corrections. In this paper, we have left this issue completely aside. However, adding real 
radiation and performing an integration over the phase space could certainly be studied in the 
effective theory framework. 

Finally, we would like to mention the possibility to extend this analysis beyond leading order 
in A2/A1. As mentioned above, if we keep only the leading terms in A2/A1 after integrating 
out the hard modes we recover the DPA of the underlying theory. Since there is nothing that 
prevents us to go beyond the leading terms in A2/A1, the effective theory approach opens up 
the possibility to go beyond the DPA in a systematic way. 

In summary, we studied the structure of higher-order corrections in the presence of unstable 
particles, exploring an effective field theory approach to the problem. Even though we made 
several simplifying restrictions this approach is very promising and seems to allow to look at a 
whole set of problems from a different, more promising, point of view. 



Acknowledgements 

VAK thanks the Leverhulme Trust for a Fellowship. Work supported in part by the EU Fifth 
Framework Programme 'Improving Human Potential', Research Training Network 'Particle 
Physics Phenomenology at High Energy Colliders', contract HPRN-CT-2000-00149. 



39 



A Expansion of loop integrals 



In this appendix we consider an example of the separation of factorizable and non-factorizable 
corrections to all orders in A2/A1. We use an expansion of loop integrals in dimensional reg- 



ularization by letting the loop momenta be in all relevant momentum regions p3| , p6[] . This 
expansion is equivalent to our effective field theory framework. 

The diagram we consider here is shown in Fig. |^(c) and contains both factorizable and 
non-factorizable contributions. The scalar integral associated with this diagram is 

dl 1 

P{P-2kl){P-2pl + D)' ^ 

where A;^ = m^, = M^, M ^ D ^ m. The integral is IR and UV finite and, thus, does not 
require regularization. The exact result is 



2 r ^2 

In — - In 



-D 



M'^-D 



M2 



TT 

"6 



2n 



(A2) 



In order to expand the loop integral in D/M"^ we should assume the loop momentum to 
be in one of the following two regions: soft, / ~ -D/M, and hard, / ~ M. Then we have to 
expand the integrand and perform the integrations over the full momentum space in dimensional 
regularization. 

' ' t(^r- (A3) 



'soft 



hard 



+ 



P {-2kl){-2pl + D) J P {P - 2kl){P - 2pl) ^^\P - 2pl 



There is only one (the leading) term in the soft momentum expansion because in higher orders 
the P term in the denominator will be cancelled. The resulting integrals are scaleless and vanish 
in dimensional regularization. The soft integral is UV divergent, and the hard integral is IR 
divergent. These divergences should be regularized by dimensional regularization and they will 
cancel in the sum. 



The soft integral is 
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For n = the hard integral is given by 
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P {P - 2kl){P - 2pl) 
whereas for n > 1 the integrals are finite and given by 
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As expected, the sum of /soft and /hard reproduce the full result. 



(A4) 



(A5) 



(A6) 
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B Matching calculation 



In this appendix we give some details concerning the matching of the {tjjtp)F^ operator defined 
in Eq. (0). To do this matching, we first have to evaluate the contribution of the operator 
(ip^)F'^, Eq. (|63D, to three-gluon processes. The Wilson coefficient 'E^^^^_^^^j^^^,^^,^(k, k'\li,l2,h) 
is then determined such that the full result in the effective theory is the same as in QCD. 

The operator {^ip)F'^ yields two different contributions to three-gluon processes. The first 
is obtained by taking the leading term in the path ordered exponentials, U{zi,Z2) = 1, while 
taking the abelian part of one of the two field tensors and the commutator term of the 
other. The evaluation of these terms is straightforward. The second contribution is obtained 
by taking the abelian part in both field tensors but the order g term in one path ordered 
exponential. In this case we have to evaluate expressions of the form 

J dxdyd'^zd'^S ^^^^i^-^) ^'^'^iiy-X) ^is3{zi-X) ^is4,{z2-X) ^-ikx ^-ik'y ^-Uizi ^-il2Z2 

S£li^2.2(^i' ^2^3, s,){-zg) r du-^ e-*- + {1 ^ 2 ^ 3} 

J zi 

where Zj, ai and are the momenta, Lorentz indices and color labels of the gluons respectively. 
The first point to note is that this contribution depends on the chosen path from zi to Z2- This 
path dependence has to be cancelled by the contribution of the operator {'il}il>)F^. From this we 
conclude that the precise form of the corresponding Wilson coefficient S*^^-* is path dependent. 

To illustrate this point we will consider the matching for two different paths. It is convenient 
to split the tensor T.^^^ into a color symmetric and antisymmetric piece 

SSll.2.2.3.3(^' ^1^1' ^2, h) = ^^;Z.U2,sJ^^ ^1^1' ^2, h) + Sgvi2.2M3.3(^' ^1^1' ^2, h) (B2) 

The tensor will be proportional to a color factor 1/2 (Tr (T^^T^-'T"-^) + Tr (T'^iT^^T"^)) 

while the S^^^") terms will come with a color factor 1/2 (Tr (T^iT^^T"^) - Tr (T'^^T'^aT^^)). 
The symmetric part is exactly the same as in the abelian theory and is given in Eq. (|6T|) with 
iV = 3. The antisymmetric part is the new feature of the non- abelian case and contains all the 
path dependence. 

The first path we consider is the linear path. In this case 

rZ2 Z^'^ Z^"^ / \ 

/ da;"^ e-^'^" = i — ^- ^ f e"*'^"^ - e^^'^^M (B3) 

Jz^ h ■ {Z2 - Zi) ^ ' 

Thus the contribution of the operator (■?/;'?/') to three-gluon processes given in Eq. ( |B1| ) will 
contain logarithms with ratios of scalar products of momenta as arguments. Since the QCD 
result does not contain such logarithms, they will have to be cancelled by the operator (■?/'?/') F^. 
As a result, the corresponding Wilson coefficient S*^^'~) will contain some logarithms. 
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In order to facilitate the presentation of the result let us introduce the shorthand notations 

= k + + (B4) 

Skk'ih, ^j, 4) = k ■ li k' ■ Ij + k ■ Ij k' ■ Ij + k ■ Ij k' ■ 4 (B6) 
and the auxiliary function /aux 

UUk, k'- h, k, h) ^ 4^.^^^.^^^^. J^^,.^^ f^,.i^ (B7) 



2 k - h k-li2 k- /i23 k' ■ h Skk'ih, h) sL'(4' ^2, ^3) 



k-h ( 1 1 \. ( k-h 

+ n . ^ ^ /, , , N log 



k ■ I123 \Skk'{h,k) slki{h,h,h) s\ki{h-,h) Skk'{h-,h-,h) ) \k ■ lu J 
The explicit form of the matching coefficient can then be written as 

^At%iM2i^2M3i^3(^' ^ 1^1' ^3) = ^ {k/iikj,^ — k^i^^ky^){k^^k^^ — k^^ku2){k^.^k^^ — k^^k^^) (B8) 

X ^/aux(^) k'] li, I2, I3) — /aux(^; k' ] li, I3, I2) — f&m{k\ k] li, I2, I3) + /aux(^'; k] li, 1^, I2) 

We should mention that the form of the matching coefficient is not unique, since it only enters 
in a symmetric combination in any calculation. The important point is that the explicit form 
of the matching coefficient Eq. (|B8|) is compatible with the estimate given in Eq. (|68|). 



Let us now consider a second path. In order to evaluate the expression given in Eq. ( Pl|) 
we choose the following path. We start at zi and follow a linear path to the intermediate point 
^oL-i _ ^as _j_ . _ zi)/l^ ■ Z3 from where we follow again a linear path to Z2- This path 
has the advantage that 



-1 '3 ■ ^3 ^3 ■ ^3 



and, therefore, Eq. ( |BT| ) will not yield logarithms upon integration. However, the presence of 
terms proportional to l'^'^ in Eq. (|B9| ) results in a more complicated tensor structure of the 
Wilson coefficient. Instead of being proportional to fc^, or k^. only, the matching coefficient will 
also contain tensor structures with one of the fermion momenta replaced by a gluon momentum. 
The explicit result is 

^^Xl>^.,sJk,k'\hj2j3)= (BIO) 
23 (k/iikl^ — k^_^ky^){k^2k'u2 ~ k'^^ky2){k^.^k'^,^ — k'^^k^^)!^^^'^ \k,k'\li,l2,l3) 

+ ^ E (-1)'^ h,X.{k,2kl2 - k',2ku2){k,X, - k'^Xs)^^'-'''^~\k,k'\k,l2,l3) 



42 



where the sum consists of twelve terms that are obtained by taking all permutations of {1, 2, 3} 
as well as the interchange terms k k'. The factor (—1)'^ denotes the sign of the permutation 
of {1, 2, 3}. The explicit form of the scalar coefficients is 

i:'-'-'-'Hk,k'\h,k,h)= (Bii) 

2k-l2 k-h k' k' -k \k-kk- k' ■ lu k ■ hj 

S(3^^^")(fc,fc'|/i,/2,/3) = (B12) 
1(2 1 \ 



3li-hk-kk'-k\k- k k' -kk' -k k- k' ■ h k' ■ I 



12 , 



Note that all terms in S^'^'"^ that are proportional to one of the gluon momenta are path 
dependent and will be cancelled by corresponding terms coming from the (■?/'?/') operator. 
Therefore, as for the linear path, the explicit form of the matching coefficient is in agreement 
with Eq. (H). 
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